Density-functional theory for ID harmonically trapped Bose-Fermi mixture 
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We present a density-functional theory for the one dimensional harmonically trapped Bose-Fermi 
mixture with repulsive contact interactions. The ground state density distribution of each com- 
ponent is obtained by solving the Kohn-Sham equations numerically based on the Local Density 
Approximation and the exact solution for the homogeneous system given by Bethe ansatz method. 
It is shown that for strong enough interaction, a considerable amount of fermions are repelled out 
of the central region of the trap, exhibiting partial phase separation of Bose and Fermi components. 
Oscillations emerge in the Bose density curves reflecting the strong correlation with Fermions. For 
infinite strong interaction, the ground state energy of the mixture and the total density are consistent 
with the scenario that all atoms in the mixture are fully fermionized. 



PACS numbers: 03.75.Mn, 71.15.Mb, 67.85. Pq 



I. INTRODUCTION 

Ultracold atomic gases provide a highly controllable 
testing ground to study fundamental problems in quan- 
tum many body physics [l| and many experimental ob- 
servations can be compared directly with exactly solv- 
able theories. The degenerate quantum gases with multi- 
component in low spatial dimensions, especially in one 
dimension (ID), become one of the growing interest- 
ing topics [2\. Multi-component gases can be mixture 
of the same species of atoms with different hyperfine 
states, i.e. spinor condensate, or mixture of different 
species of atoms. The competition between the intcr- 
and intra-species interaction makes the mixture system 
more complicated and there exhibit richer physical phe- 
nomena than its single-component counterpart. Bose- 
Fermi mixture is originally realized in experiments as 
a result of sympathetic cooling technique, i.e., cooling 
the fermions to quantum degeneracy through the me- 
diation of bosons tSnZi]- Then many theoretical studies 
have been performed on three dimensional mixture, deal- 
ing with the phase separation [a|, pairing Q, superfluids 
and Mott insulators transition [l3|, BEC and Bardeen- 
Cooper-Schrieffer crossover [il|, etc. On the other hand, 
ID systems attract attentions for the simplicity of the- 
oretical models and for the significance of quantum cor- 
relation effects therein [12]. Experimentally the ID sys- 
tems can be realized by confining the cold atoms in two 
dimensional optical lattices or in strong anisotropic mag- 
netic trap [i3l ■ The interaction among the atoms can be 
tuned in the whole regime of interaction strength via the 
magnetic Feshbach resonance and controlling the trans- 
verse confinement of magnetic trap [l4|. Interestingly 
enough, the properties of these system, including the 
ground state, elementary excitations as well as thermody- 
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namics, are sometimes fairly well captured by the exactly 
solvable models studied a few decades ago [l^l ■ 

Many theories have studied the ID Bose-Fermi mix- 
ture both homogeneous and trapped gas in external po- 
tentials. When the system is homogeneous. Das [Ig] early 
plots the ground state phase diagram based on the mean- 
field theory, which predicts the occurrence of phase sepa- 
ration (i.e., demixing) of the two components. Luttinger 
liquid formalism shows that for strong enough repulsion 
the two components of the mixture, with bosons either 
a quasicondensate or impenetrable particles, repel each 
other sufficiently to demix J^]. But the exact Bethe 
ansatz solution for the ID mixture with equal mass and 
equal coupling constants points out the absence of demix- 
ing [1^ [I^ . It indicates that mean-field theory and Lut- 
tinger liquid theory are reliable only for very weak inter- 
action. With the system loaded in optical lattices, the 
phase diagram and correlation functions have been in- 
vestigated with the bosonization method 20] and quan- 
tum Monte Carlo (MC) simulation 2\^- For the trapped 
system, the local density approximation (LDA) on the 
Bethe ansatz solution shows that the harmonic trap ID 
mixture would partially demix for strong repulsive inter- 
action [l9|. The finite temperature Yang- Yang thermo- 
dynamics and the quantum criticality analysis based on 
thermodynamic Bethe ansatz (TBA) both support the 
description of phase separation [22]. In the infinitely 
strong interaction limit, i.e. the Tonks- Girardeau (TG) 
regime, Bose-Fermi mapping method [23| gives a result 
that the density profiles display no demixing among the 
two component, where the exact ground state is hi ghly 
deg enerate and the most symmetrical one is chosen [2J- 
l2a ]. Later detailed calculations are done for all degener- 
ate manifolds of the ground state [27] and the conclusion 
of nondemixing remains for the mixture in TG limit. 

So far a method is not available for ID trapped Bose- 
Fermi mixture suitable for the whole repulsive inter- 
action regime. Numerical simulations such as Density 
Matrix Renormalization Group (DMRG) and quantum 



MC are however limited to few atom numbers in lattice 
models [2l|, [23|- In this paper, we develop Hohenberg- 
Kohn-Sham Density-functional theory (DFT) to inves- 
tigate the ground state properties of ID harmonically 
trapped Bose-Fcrmi mixture. It is well known that DFT 
is a successful and widely used approach for treating 
the electron systems with long range Columb interaction 
[28, 29] . Recently it has been successfully generalized to 
cold atom systems with short-range contact interaction 
for three dimensional bosonic atoms [30, [31| and three di- 
mensional Bose- Fermi mixture [32| • In the framework of 
DFT, in order to investigate the ground state properties 
of an inhomogeneous interacting system, a homogeneous 
interacting system is often needed in the process of local 
density approximation (LDA) for the exchange correla- 
tion energy [2S]- Because Bethe ansatz method can give 
exact solution for ID homogeneous systems, several au- 
thors developed the DFT based on Bethe ansatz results 

to solve the ID Bosons ^U, ^3, H| and ID Fermi cold 

I — I \ h — — — 
atom systems [35|, 136| . Here we apply this method for 

ID Bose-Fermion mixture. As can be seen below, the key 

points of our scheme include a suitable fitting formula for 

the ground state energy and appropriate choices of the 

functional orbitals for boson and fermions. 

The paper is organized as follows. In Sec. II we derive 

the universal Kohn-Sham equations with LDA for ID 

Bose-Fermi mixture and then present the expression of 

exchange-correlation energy. From the exact result of 

Bethe ansatz for homogeneous gas with equal masses of 

atoms and equal interaction of boson-bosn and boson- 

fermion we find a fitting formula for the ground state 

energy to simplify the numerical iterations. Then the 

equations are solved numerically and the ground state 

energy and density distribution are discussed in Sec. III. 

Finally we conclude our result in the last section. 



Here tub, mp are boson and fermion masses, Vb{x), 
Vf{x) are external potentials, and gsB, gsF are the ef- 
fective ID Bose-Bose and Bose-Fermi interaction param- 
eters, which can be tuned experimentally [ijj 123 • The 
Fermi-Fermi interaction is not considered because the 
Pauli exclusion principle suppresses the contact s-wave 
scattering and their p-wave scattering can be neglected. 

According to the Hohenberg-Kohn theorem I of DFT 
[28| , the ground state density of a bound system of in- 
teracting particles in some external potential determines 
this potential uniquely. It thus gives us the full Hamilto- 
nian ([T]) and particle number TV. Hence the density deter- 
mines implicitly all properties derivable from H through 
the solution of the time-independent or time-dependent 
Schrodinger equation. Though proved originally for 
Fermions, the theorem can be straightforwardly gener- 
alized to Bosons as well as the mixture of Bosons and 
Fermions studied here. Denote the densities of bosons 
and fermions as ub {x) and np {x) , respectively, the to- 
tal density is then obviously n (x) — ub (x) + np (x) . 
The number of Bosons and Fermions are conserved sep- 
arately, i.e. JnB{x)dx — Nb, J np (x) dx — Np and 
J n (x) dx = N. The ground state energy, defined as 
{g\ H \g) with \g) the ground state of system, is a func- 
tional of the densities Eq [ub {x) , np (x)] , which can be 
decomposed as 
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II. THEORY 

A. Kohn-Sham equations 

We consider a ID trapped mixture of Nb bosons and 
Np spin-polarized fermions with two-body contact inter- 
actions. N = Nb + Np is the total atom number. The 
system is described by the Hamiltonian 
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The first two terms are Bose and Fermi kinetic energies of 
a reference noninteracting system. The next two terms in 
the second row are external potential energies and those 
in the third row are Hartree-Fock energies (i.e., the mean- 
field approximation of the interaction energy). The last 
term is the exchange correlation energy which includes 
all the contributions to the interaction energy beyond 
mean-field theory. 

We assume the Bosons are in a quasi-condensate state 
and Fermions are in normal state. Thus we introduce a 
single Bose functional orbital (p (x) and Np Fermi func- 
tional orbitals ipj (x) {j ~ 1 ■ ■ ■ Np) which are orthogo- 
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This is different from the way of 
where only one condensed orbital of fermionic 
pair is considered for the mixture of bosons and paired 
two-component fermions in superfluid state or BCS state. 
With (l>{x) and ipj (x), the densities are expressed as 

nB (x) = Nb4>* {x)(t){x) , 

Nf 

np (x) = Yl ^J (^) ^J (^) ■ (^) 



and the kinetic energies are 
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As far as the exchange correlation energy E^c [ns, np] is 
concerned, when the confinement is weak, we adopt the 
Local Density Approximation (LDA), i.e., the system can 
be assumed locally homogeneous at each point x in the 
external trap. In this way Exc is approximated with an 
integral over the exchange-correlation energy per atom of 
a homogeneous interacting mixture £\°™ {riB, np) 
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where the densities hb, np are taken at point x. 
For such a homogeneous interacting mixture. 
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where ehom jg ^^^^ ground state energy per atom; s\°™ — 
gBBn%/2n + gBF^BnE /n is the mean field interaction 



energy per atom; k^°™ = ft^TT^n^/GmFn is the kinetic 
pressure terms, i.e., the total kinetic energy dividing by 
the total number of fermions and bosons in a noninter- 
acting homogeneous mixture. Here the kinetic energy 
of the bosons is easily shown to be zero and the kinetic 
energy comes solely from the exclusive quantum state 
occupation of fermions . 

Hohenberg-Kohn theorem II [2^ guarantees that the 
ground state density distributions is determined by varia- 
tionally minimizing Eq with respect to hb (x) and np (x) , 
which is equivalent to a variational calculation of ^ with 
respect to the Bose and Fermi functional orbitals 0* , ip* 
After substituting (O into dH) and substituting (jH])- 
into ([2]), we carry out the functional derivatives 



SiEo- eNB{ / (j)*dx(j} - 1) 1 / 50* = 0, 
sIeo-J^Vj {j r,dx^, - 1) / 5^* = 0, (8) 



where e and rjj {j ~ 1^2,- ■ ■ Np) are Lagrange multipliers 
conserving the normalization of 4>{x) and ^j (x). Then 
we can get the Kohn-Sham equations (KSEs) 
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Here a^b "" = d (ne*^""") /Oub and 
d (n£^°™) /dnp are Bose and Fermi chemical potentials 
of a homogeneous interacting mixture. Physically e and 
11 j are the lowest eigenvalues of KSE. In ([3|), the sum in 
np (x) runs over the occupied orbitals ipj with lowest rjj. 
Left multiplying ip* on both sides of pop . performing 
summation over j and integrating over x, we get an ex- 
pression of Tp'^-' defined in ([5]). Analogously from the 
normalization of 4>{x) we may get an expression for T^'^ 
defined in ([4]). Inserting these two kinetic terms into Q, 
the ground state energy ^ is expressed as a function of 
e and rjj 
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If e"°™ [nB,np] are known, we can solve the KSEs to- 



gether with ^ to find the density distributions nB (x), 
np (x) and then calculate the ground state energy Eq 
from (jlip. In the following we present the result of 
^hom ^^g^^pj \)y means of the Bcthe ansatz method. 

B. Ground state energy of homogeneous system 

In the absence of an external trap the system is homo- 
geneous, which can be solved exactly via Bethe ansatz 
method for a much restrictive but simple case 
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It describes the situation that the interactions of Boson- 
Boson and Boson-Fermion are repulsive with the same 
strength, and the masses of Boson and Fermion are the 
same too. Detailed possible ways to realize this situa- 
tion in cold atom experiments have been considered pre- 
viously |19j . The first condition can be satisfied using 
the combination of Feshbach resonance (to control the 



interactions) and appropriate choice of the tuning of the 
trapping laser frequencies (to adjust the the ratio of the 
radial confinenient of Bosons and Fermions). The sec- 
ond condition is approximately satisfied with a mixture 
of two isotopes of a species of atoms. Isotope mixture 
is widely used in experiment for it can avoid the gravi- 
tational sag of an external potential caused by different 
masses. The experiments have realized 3D isotope mix- 
tures ^Li-^Li 0, i^^Yb-i^^Yb d, 40k;_41j^ Q ^^^^ ^g g^g 

no obvious obstacles in ID. Under these two conditions, 
the Hamiltonian of ID homogeneous Bose-Fermi mixture 
is 
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This model is solved by means of Bethe ansatz method 
by Lai and Yang in 1971 for the ID mixture of bosons 
and spin-1/2 fermions 18]. Imambekov and Demler in- 
vestigated the ground state properties in detail for the ID 
mixture of bosons and spin-polarized fermions [l9[ and 
extensive studies have been done in [2^, [39l44l| includ- 
ing the thermodynamics and correlation functions. Here 
we briefly review the main results of [1^ [l^ |4l[ which 
are readily used as the homogeneous reference system in 
our DFT theory. Under the periodic boundary condition 
and in the thermodynamical limit (the system size and 
the number of atoms are infinitely large but the atomic 
densities are kept finite), the ground state Bethe ansatz 
integral equations are 
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where c = mg/fi^, k and A are the quasi-momenta and 
spectral parameters and p (k) and a (A) denote their cor- 
responding density distributions. The integration limits 
B and Q are determined by the normalization condition 
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The ground state energy per atom is written in our no- 
tation as 
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For convenience, let us define the fraction of Bosons 
a = nBJn and the dimensionless Lieb-Liniger parameter 
7 = mgj {Y?n) . We then introduce variables x — k/Q 
and y = kjB such that p (fc) = p ixQ) = gc ix) and 
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9c [x) = — 

ZTT 



gs (y), and (|14 p - (fT6|) are transformed 
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Here the function 



e(7,a) = ^ 



x^gc (x) dx 



(18) 
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can be solved numerically with the combination of (|17p 
and (|18l) by the iteration method. In the limiting cases of 
a = 0, 1, the system is purely fermions or purely bosons, 
e (7, 0) — 77^/3 is a constant while e (7, 1) coincides with 
gl-l (7) in the Lieb-Liniger model [42|. When the in- 
teraction is weak, 7 *C L, the mean field result of p^ 
is already available in [la, [l3|; when the interaction is 
strong, 7 :^ 1, one can neglect the dependence of the 
first integrand in ([T7|) on x and gc (x) can be approxi- 
mated as a constant gc- Therefore we get the asymptotic 
behavior of the function e (7, a) for 7 
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e (7 — > +co,a) — 
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where F (a) = a + sin {air) /n. In the limiting case of 
0, e (0, a) = TT^ (1 - a)^ /3, therefore e^""" (n, 0, a) 
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the energy comes solely from the 



kinetic energy of free fermions. In the Tonks-Girardeau 
limit, e (-l-c», a) = tt^/S, which means the energy of Bose- 
Fermi mixture with infinitely strong repulsive interac- 
tions is equal to the energy of all atoms treated as free 
fermions. 

For practical use, we need calculate e (7, a) for a lot 
of points (7, a). If we use numeric iteration method for 
every point, it will be very time-consuming. To avoid 
this, we managed to retrieve a parametrization formula 
for e (7, a) based on the above limiting cases, which reads 
as 
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-(l + /2(7))(l-a)'). 
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Here /i (7) is the approximation of g^-l (7)- We give 
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(23) 



with 61 = 11.37, 62 = 4.68, ai = 12+6i-462, a2 = -4+^2 
and 60 = 7r^ai/3, which exhibits the same asymptotic 
behavior as e^-^ (7) in the weak and strong interaction 
cases to the order of 7 and 1/7^, respectively. The func- 
tion /2 (7) is determined by the numerical iteration result 
for some sampled values of 7, and we fit it as 

/2 (7) = Co exp(ci7) - (co + 1) exp(c27) (24) 

with Co = 0.21, ci — —0.02, C2 = —1.45. e (7, a) gives 
the exact behavior at the limits a = 0,7 = 0,+oo and 
approximates to bl-l (7) at the limit a = 1. In in- 
termediate value of a and 7, e (7, a) deviates with a 
maximum relative error of 0.03 from the numerical re- 
sult at 7 sa 2.5, a ss 0.9. In Fig[Tl we exhibit the result 
of exact numerical result of e (7, a) compared with the 
fitting formulas e (7, a) for various interaction strength 
7 and the fraction of Bosons a. Clearly the fitting for- 
mulas represent quite well the Betha-Ansatz result for 
the whole range of interaction and arbitrary fraction of 
bosonic atoms in the mixture. These formulas are then 
adopted in the following solution of the KSEs equations. 
From ehom (^j^^q,^^)^ the ground state Bose and Fermi 
chemical potentials can be obtained as 
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where 



fF{l,a) = 3e-7— -a— , 



(26) 



with /B(0,a) = 0,/F(0,a) = tt^ (1 - a)^ and 

JB (+00, a) = fp (+00, a) = TT^. 



III. NUMERICAL RESULTS 
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FIG. 1: (Color Online) The function e (7, a) for the ground 
state energy of homogeneous Bose- Fermi mixture system. Nu- 
merical exact result (red dashed line), obtained from the so- 
lution of p7|l and (|18|l is compared with the fitting formula 
(solid black line) given by ^. (a) 7=0,0.5,2,5,10,20,100 from 
bottom to top; (b) a=0, 0.05, 0.1,0. 2, 0.3,0. 5,1 from top to bot- 
tom. 



DFT result, we first discuss the KSEs for some limiting 
cases. 

When there is no interactions in the mixture, [/ = 0, 
KSEs correctly reduce to the equations for noninteract- 
ing Bosons and noninteracting Fermions in the harmonic 
trap. The densities of bosonic and fermionic components 
are respectively 



Inserting (gS]) into the KSEs ^ and (HU]), and assum- 
ing Bosons and Fermions suffer from the same harmonic 
external potentials Vb{x) — Vp{x) = muj'^x^ /2, with lo 
is frequency, we can get the ground state density profiles 
of each component by solving the KSEs together with 
the constraint ([3]) by means of numerical iteration. The 
ground state energy follows immediately from ([Tl]) . Here 
we introduce the length unit a = ^Jhjrmjj and a dimen- 
sionless interacting parameter U = g/ahu such that the 
space dependent Lieb-Liniger parameter is expressed as 
7(2;) = U/an{x). Before going into the details of the 
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where Hi (x) is the Hermite polynomials. Here the nonin- 
teracting Bose density profile (see the topmost black line 
in Fig. [2]) is a standard Gaussian-like shape and Fermi 
density profile (see the red dotted line in Fig. ^ is char- 
acterized by a half ellipse-like shape with Np oscillations. 
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FIG. 2: (Color Online) Density distribution of Nf — 10 non- 
interacting Fermions (red dotted line) and the density distri- 
butions of Nb = 10 Bosons for different interaction parameter 
U. The six black solid lines from top to bottom are respec- 
tively for Bosons with f/ = 0, 0.1, 0.5, 1, 5, 10, +00. 



50 
40 
30 
20 



3 








-10 

-20 
-30 



/ 
/ 


y 


^ 


-- 




// 


/ 






■ 


/ 


X" 




— 


" 




1' / 










/ 








■ 


\ 








F .... T^'f 






tg ■■ ■■ 1 










■■E^^ — -i?"' ■ 






\ 








—-Ef" ' 



1 
16 


/ 


14 

12 

10 

8 


/ / 
/ / 

/ ' 
/ / 
/ / 
/ / 
/ / 


6 


/ / /" 


4 




2 


1 - ■ . 





... 


-2 


'■•• 



10 



20 30 
U 



40 



50 0.0 0.5 1.0 

u 



FIG. 3: (Color Online) Evolution of energies of Nb = 10 
Bosons with increasing interaction parameter U . Contribu- 
tions to the ground state energy Eo include: kinetic en- 
ergy T*"*^', external potential energy Epot, Hatree-Fock energy 
Ehf, exchange correlation energy E^c- Right Panel: Details 
in the mean-field regime. 



The ground state energies of these two components are 

(29) 

EoF - Y. V+2l^ ^^°^ 
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and the total ground state energy is Eq = Eqb + -Eqf • 

When the interaction is weak, neglecting E^c in ^, 
the KSEs reduce to our familiar mean-field formulas 
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For bosons it is nothing but the Gross-Pitaevskii equation 
for dilute gas. The equation for fermions, on the other 
hand, reminds us the superfluid theory of the mixture 
of bosons and paired BCS states of the two-component 
fermions where only one Fermionic orbital is considered 

When the interaction is strong, for system of large 
atom numbers Nb,Nf ^ 1, one can safely use the 
Thomas- Fermi approximation (TFA), i.e., the kinetic en- 
ergies r^^ and r^^ in the energy functional ([2]) are ap- 
proximated to zero and J n (x) Kg°™ (x) dx, respectively. 
Minimizing Eq directly with respect to ub (x) and np (x), 
we get the TFA formulas 



mwV-KAi^°"^([nB,n^];x) = /.t°,, (32) 



where /i^ and /i^ are constants fixed by the normaliza- 
tion conditions / hb (x) dx = Nb and J np (x) dx ~ Np. 
Eqs. (j31]) are explained as the LDA of the chemical po- 
tentials at point x in [19j and have been used extensively 
[22, jlOj. That means, in slowly varying external har- 
monic trap chemical potentials at point x are related to 
those in the trap center a; = (/x^ and fi^p). 

When repulsive interactions are infinitely strong, 
^hom ^ ^hom ^ h'^n^n^ /2m, ^ reduces to a single 
equation 
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with /i'^ is decided by J n (x) dx = N. This gives us the 
explicit result for total density distribution 



n (x) = 
and ground state energy 
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We see that they are exactly the density distribution and 
energy of N free fermions in a harmonic trap. Equation 
([33| . however, gives nothing about the densities of Bose 
and Fermi components. The method here is insufficient 
for the infinitely strong interaction. We may, on the other 
hand, resort to the Bose- Fermi mapping method [23, 123] 
which gives the Bose and Fermi density profiles as 



/ \ ^B,F , 2/2 

nB,F [x) = -rn= cxp (-X a 



E 

ri=0 



2"n! 



(36) 



The two components are nondemixing in agreement with 
the generalized Bcthe ansatz wave function |27| . 

The DFT results are summarized in Figs. 2-8. First, 
for a pure bosonic system, equation ^ is just the general- 
ized Gross-Pitaevskii equation appeared in Refs. [3ll.l43|. 
We show the density profiles for Nb = 10 bosons in Fig. 
[5] for the cases of [/ = 0, 0.1, 0.5, 1, 5, 10 and +oo, respec- 
tively. With the increasing of U, the density profiles vary 
from a standard Gaussian-like to a non-oscillating half- 
ellipse shape. Comparing the density profiles of Bosons 
at U = -|-(X) and the noninteracting fermions, we find 
that they match each other quite well except the density 
oscillations. The results mean the density distribution of 
Bosons with infinitely strong repulsive delta interaction 
is basically the same as that of a noninteracting Fermi 
gas which is consistent with the theory of Bose-Fermi 
mapping theory [44j . Our theory fails to reproduce the 
density oscillation due to the impenetrable property of 
ID system with strong interaction because we adopt one 
single functional orbital (p{x) for the ID Bose liquid. The 
exact oscillations reflecting the structure of the occupied 
orbitals should be quested from the real wave function 
such as by the exact diagonalization method [45| . In the 
limit of large particle number the differences between the 
oscillating and non-oscillating profiles become impercep- 
tible. 

The ground state energy evolution as a function of U 
is illustrated in Fig. [31 We can see that with the increase 
of U, the kinetic energy T^'^f decreases slowly indicating 
that the interaction restrains the movement of atoms. 
The external potential energy Epot increases as a result 
of wider and wider occupied regime of the trap. Both 
of these two energies evolve to constant energies. The 
Hatree-Fock energy Ehf increases almost linearly while 
the exchange correlation energy E^c decreases in the 
whole interaction regime. These two terms play more and 
more important roles in the DFT theory for stronger in- 
teraction and they approximately cancel each other. All 
these energies contribute to the total energy Eq, which 
starting from the noninteracting value 5huj approaches 
the strongly interacting limit 50haj. For U < 0.9, the 
exchange correlation energy is much less than the total 
energy, \Exc/Eo\ < 0.1, which can be seen as the effec- 
tive regime of mean field theory. For a TG gas with 
U — -l-oo and a chemical potential /i^°™ = h^ii^n? /2ra, 
numerically solving equation ^ gives Eq = 50.5024^10; 
which lies slightly above 50fia; because the introduced 
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FIG. 4: (Color Online) Density distributions of a mixture of 
Nb ~ 10 bosons and Nf — 10 fermions for different inter- 
action parameter U = 0.1,1,5 (left) and U — 20, 100, -l-oo 
(right). 




X [a] 

FIG. 5: (Color Online) The total density profiles n as a func- 
tion of X for a Nb = Nf = 10 mixture with U = -|-cx). The 
result of KSEs (black solid line), TFA (blue dashed line) and 
Bose-Fermi Mapping (red dotted line) are compared and the 
inset shows a zoom into the structure of the ocsillations. 



Bose functional orbital ^(a;) is only an assistant varia- 
tional function instead of the true wave function of the 
interacting Bose system. 

We now turn to illustrate the main result of a mixture 
of Nb = 10 bosons and Np = 10 fermions. The densities 
of non-interacting mixture ([?7)) and (|28p are taken as the 
starting point of the iteration of KSEs (|9ll0p for a small 
interaction parameter (e.g. U = 0.1). The eigenvalues 




FIG. 6: (Color Online) Ground state energies as a function of 
U for mixture of Nb = Nf — 10. Contributions to the ground 
state energy Eo are similar to those in Fig. 3. All terms but 
the exchange correlation energy originate from bosons and 
fermions. Right Panel: Details in the mean-field regime. 



e, rjj and functional orbitals 0, ipj are found by iterating 
to the desired degree of accuracy. The new densities are 
initial densities for the next iteration for a larger interac- 
tion parameter, and so on. The density profiles for differ- 
ent U are displayed in Fig. S] It shows that with increas- 
ing U, the peak of the total density n{x) decreases niono- 
tonically and atoms tend to occupy wider regime. The 
density of Fermi component changes smoothly in ampli- 
tude, while Bose component becomes more and more flat 
and ripples begin to appear for stronger interaction. At 
weak interaction, C/ = 0.1, 1 as in Fig. |4](a) and (b), both 
Bosons and Fermions are located in the center of the trap. 
For an intermediate interaction strength, e.g. L'^ = 5 as in 
Fig. 111(c), some fermions are excluded from the center of 
the trap while Bosons are held mainly in the center. We 
notice that oscillations emerge in the Bose density curves 
reflecting the strong correlation with Fermions. When U 
becomes further stronger, U = 20, 100 as in Fig. Illjd) 
and (e), more Fermions are repelled out from the center 
and a clear signature of phase separation of bosons and 
fermions is seen in the figures. High density of discrete 
bosons are surrounded by fermions, which nevertheless 
still have chance to squeeze between the opening space of 
bosons. The total density profile approaches half ellipse- 
like for U = -t-oo as shown in Fig. IHf). We may have 
a close inspection of the case of infinitely strong inter- 
action. In Fig. [5l the DFT result of the total density 
is compared with those from TFA and Bose-Fermi map- 
ping. It is clear that the agreement is fairly good except 
tiny difference in the number, position and amplitude of 
the oscillations which are enlarged in the inset of Fig. \5\ 
Again for large atom number the differences between the 
oscillating and non-oscillating curves is unperceivable. 
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FIG. 7: (Color Online) Density distributions of mixture with 
A^ = 20 for different mean fraction of bosons a — Nb/N and 
different interaction parameter U. 



Fig. ini describes the evolution of all contributed en- 
ergy terms in ([2]) as a function of U. The line types 
denote the same energy terms as in Fig. [3] except that 
here the kinetic energy and external potential energy re- 
spectively include two terms relating to Bosons (in red) 
and Fermions (in blue). The trend of these lines re- 
semble those in Fig. [3] for the same reason. Especially 
the exchange correlation energy here contest with two 
Hatree-Fock terms representing the mean field energy of 
boson-boson and boson-fermion interaction respectively. 
The total ground state energy evolves from 55huj, the en- 
ergy of 10 ideal bosons and 10 ideal fermions, to 200haj, 
the energy of 20 fully fermionized atoms according to 
(l29t and ([30]). In the parameter range of < [/ < 2, 
\Exc/Eo\ < 0.1, mean field theory is regarded effective. 
At U — +CXD, we numerically obtain an upper bound 
for the ground state energy Eq — 200.1364?ki; which 
is very close to the exact result of full fermionization 
Eo = 200fiaj. 

There are some two contradicting predictions for the 
spatial structure of the components densities of trapped 
TG mixture. Based on the exact result of Bethe ansatz 
and TFA formulas ^ (where they called it LDA), Ref. 
(l9| gives a phase separation result of the two compo- 
nents for a strong but finite interaction. The thermody- 
namical Bethe ansatz (TBA) at finite temperature [23| 
tends to support this scenario. The Bose-Fermi map- 
ping methods in Ref. [2j] and [2J] is proposed in TG 
limit of infinitely strong interaction and the results show 
that the component distributions of the Nb — Np mix- 
ture are completely the same according to (1361) therefore 
display no demixing. But the authors of [2J] have no- 
ticed that the ground state given in this way is highly 
degenerate. They subsequently used a generalized Betha 
ansatz wave function, in which the 'orbital' part of the 




FIG. 8; Ground state energy Eo as a function of U for different 
fraction of bosons a in the mixture with total atom number 



wave function is essentially replaced by a Slater determi- 
nant of single-particle Schrodinger equation in the trap 
potential, to give a nondemixing result at finite large in- 
teractions. They further tested this nondemixing result 
with numerical DMRG simulations for a lattice model of 
Nb = Np = 2 mixture [23|. We observe, however, obvi- 
ous signature of phase separation in Fig. 4 of [27| for rel- 
atively large interaction U = 100. The intrinsic nature of 
the phase separation and nondemixing in TG limit orig- 
inates from the Bethe ansatz and Bose-Fermi mapping 
techniques respectively. We expect experimental verifica- 
tion of the nature of spatial configuration about trapped 
ultracold atomic mixtures. 

Finally we discuss the effect of another system param- 
eter, i.e. mean value of the fraction of bosonic atoms 
number a — Nb/N, on the density profiles. Fig. [7] 
shows the density of each component and the total den- 
sity of iV = 20 atoms in the mixture with a = 0.25, 0.75 
and interaction parameters U — 1,10. Fig. |S] compares 
their energies. Bosons will dominate the total density 
profile when more bosons are put into the mixture for 
weak as well as strong interaction. When more fermions 
are prepared in the gas, bimodal distribution is clearly 
seen the total density where bosonic Gaussian shape is 
superimposed onto the fermionic shell-like structure. In 



the strong interaction limit the total density approaches 
to the typical half-ellipse no matter how many bosons 
or fermions are involved in the mixture. The number 
of fermions contribute to the ground state energy more 
effectively in the weak interaction case. This situation 
changes for strong interaction where the energies for all 
values of a approximate to the limit value of the fully 
fermionization of the system. 

IV. CONCLUSION 

In conclusion, using the DFT we study the ground 
state energy and density distribution of the Bose-Fermi 
mixture in a quasi ID harmonic trap. Based on the Bethe 
ansatz solution for the mixture, we managed to obtain a 
fitting formula for the function e(7, a) for the ground 
state energy of homogeneous system. The KSEs are ob- 
tained from the variational minimization of the energy 
functional of trapped mixture with respect to the densi- 
ties of bose and fermi components. We found that when 
the interaction between the atoms varies from zero to 
positive infinitely, the ground state energy of the mix- 
ture would evolve to the constants of the noninteract- 
ing Fermions and the total density approached a half el- 
lipse profile. More and more fermions are repelled out 
of the trap center, while bosons occupy the central re- 
gion. Phase separation of boson and fermion components 
occurs for strong interaction in agreement with the the 
result Bethe ansatz method plus LDA. The calculation 
here applies equally to the pure bosonic case, different 
fraction of bosons, as well as in the TG limit. Our DFT 
theory is also suitable for mixtures in optical lattice and 
could be extended to study the dynamical and thermo- 
dynamic phenomena of the mixture. 
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